Geometry and topological order in the non-relativistic Luttinger liquid. 
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We demonstrate that the notion of spin-charge separation as follows from the Tomonaga-Luttinger 
theory of relativistic fermions does not suffice to completely specify the Luttinger liquid state as- 
sociated with the fermions of condensed matter physics. The latter carries in addition a form of 
topological order which can be made visible using non-local correlation functions familiar from lattice 
gauge theory. 
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Systems composed of an infinite number of interacting 
quantum particles can be characterized by long wave- 
length excitations carrying quantum numbers unrelated 
to those of the constituents. A classic example is the 
phenomenon of spin-charge separation which is well es- 
tablished in one dimensional electron systems The 
elementary excitations can be viewed as pieces of the 
original electron, carrying separately the charge and the 
spin. This notion first emerged when it was demonstrated 
that systems of relativistic fermions in l-flD can be rep- 
resented exactly under a variety circumstances in terms 
of Bose fields, carrying separately the spin and the charge 
degrees of freedom (bosonization) . Although a- priori it 
is unclear whether the relativistic theory has anything to 
say about systems of non-relativistic electrons, an over- 
whelming case emerged that bosonization is correctly de- 
scribing the asymptotic properties of the local correla- 
tion functions of the latter. However, the loophole is 
that these might differ by a non-critical topological order 
which cannot be detected using two-point correlators in 
terms of local operators. Here we will demonstrate that 
at least the Hubbard model in 1+lD carries such an or- 
der which is alien to relativistic fermions. We believe that 
the Hubbard model is in this regard fully representative 
for all non-relativistic electron systems. 

This hidden order has to do with the fact that the 
charge carriers ('holons') are at the same time Ising 'do- 
main walls' (kinks) in the spin system, a property shared 
with stripes in two dimensional systems Resting on 
the geometrical definition of spin-charge separation as de- 
duced from the Bethe-Ansatz solution of the Hubbard 
model at infinite U , obtained first by Woynarovich Q and 
later Ogata and Shiba we construct a non-local op- 
erator measuring directly the topological order. Using a 
state of the art DMRG method, we demonstrate numer- 
ically that this order is present in the Hubbard chain for 
all positive C/'s and electron densities. We subsequently 
show that a straightforward application of bosonization 
leads to inconsistencies, demonstrating that this topolog- 
ical order is not present in the Dirac vacuum in 1-l-lD. 
Our non-local correlator has the structure of the Wil- 



son loop of a Z2 gauge theory and we conclude with the 
speculation that a specific kind of 1+lD superconductor 
might be governed by a local Ising symmetry. 

Some time ago , a simple but most peculiar prop- 
erty of the Bethe-Ansatz solution of the ID Hubbard 
model in the limit that J7 — s- 00 was found; namely that 
the wave function factorizes into a charge part ipSFi de- 
pending only on where the electrons are, and a spin part 
TpH which merely depends on the way the spins are dis- 
tributed, 

{%}fil') = ^SF[{x.}tl) V'i/({2/,}fi?). (1) 

The charge part ipsF is the wave- function of a non- 
interacting spinless-fermion system where the coordi- 
nates Xi refer to the positions of the electrons/spinless 
fermions. The spin part ij^H is identical to the wave- 
function of a chain of Heisenberg spins interacting via 
a nearest-neighbor antiferromagnetic exchange. In ipH 
only the positions of the up spins are needed and these 
correspond with the coordinates yi. However, the mira- 
cle is that the coordinates yi do not refer to the Hubbard 
chain, but instead to a new space: a lattice with sites 
at coordinates xi, 0:2, xat given by the positions of the 
charges in a configuration with amplitude ipsF- 

This shows that the quantum dynamics of interacting 
electrons generates a geometrical structure analogous to 
the fabric of general relativity. Let us visualize this for 
a representative example (Fig. 1). Consider N electrons 
on a chain with L sites under the condition that N < L 
such that the charge configurations can be specified by 
the locations of the holes. A charge configuration in the 
full Hubbard chain ('external space'), has an amplitude 
tpSF in the wavcfunction with the coordinates of the dots 
corresponding with the x^'s. The spin system sees a dif- 
ferent 'internal space' obtained from the full space by re- 
moving the holes together with the sites where the holes 
are located, substituting the hole and its site with an an- 
tiferromagnetic exchange between the sites neighboring 
the hole (the 'squeezed space' ||]). 
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FIG. 1. The emergence of the Z-i 'sublattice parity field' 
in the geometrical 'squeezing' operation. The open circles re- 
fer to some configuration of holes. The black dots define the 
embedding space of the spin system in both the full- (upper) 
and squeezed (lower) lattice. The unsqueezing operation can 
be parametrized in terms of binding of the electric charge 
quantum to flips in the Z2 valued sublattice parity field. 

This is a very simple geometrical structure which can 
be formulated in terms of a simple topological 'gauge' 
theory. In what regards are the full chain and the 
squeezed chain different? The squeezed chain is obvi- 
ously shorter than the full chain and this is a simple 
dilation: a distance x measured in the full chain becomes 
a distance px in the squeezed chain (p — N / the elec- 
tron density) when a; 3> 1, the lattice constant. The 
other aspect is also simple, but less trivial in its conse- 
quences. The spin system is a quantum-antiferromagnet 
and is therefore sensitive to the geometrical property of 
bipartiteness. A lattice is called bipartite if it can be 
subdivided in two sublattices A and B, such that all sites 
on the A sublattice are neighbored by B sublattice sites 
and vice versa. This division can be done in two ways 

( A^ B - A- B--- and B - A - B - A---) 

defining a Z2 valued quantity 'sublattice parity', p — ±1. 
Consider now what happens with sublattice parity in the 
squeezing. For the Heisenberg spin chain a redefinition 
of p = 1 ^ — 1 does not carry any consequence ('pure 
gauge'). However, sublattice parity becomes important 
in the mapping of squeezed space into full space (Fig. 1). 
'Fix the gauge' in squeezed space by choosing a particu- 
lar sublattice parity, and consider what happens when it 
is unsqueezed. The holes are inserted, and because every 
hole is attached to one site, every time a hole is passed 
the sublattice parity flips. This is true for every instan- 
taneous charge configuration, but the ground state is a 
a superposition of many of these configurations: charge 
is fluctuating and since it is attached to the sublattice 
parity flips, the full space which is observable by exter- 
nal observers (experimentalists) should be considered as 
a fluctuating geometry. However, this is a very simple 
fluctuating geometry because all that is fluctuating is the 
property of bipartiteness. 

Let us now ask the question if a correlation function 
can be defined acting on the full Hubbard chain which 
can measure the true spin correlations living in squeezed 
space. Since all that matters is sublattice parity this can 
be achieved by simply multiplying the spin operator by 
a factor —1 every time a hole is passed, thereby remov- 



ing the sublattice parity flips from the spin correlations. 
Define staggered magnetization as M{xi) = {—lY^ S{xi) 
where S is the spin operator [S^ = ^(n-f — n^), — cjcj^) 
with the charge operator n^.. = n^^^ + f^xii taking the 
values 0, 1 and 2 for an empty-, singly- and doubly occu- 
pied sites respectively. The correlation function we are 
looking for is P], 
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The operator (— l)(i^"3) takes the value -f 1 for a singly 
occupied site while it is —1 for a charge (hole, or doubly 
occupied) site. By multiplying these values on the inter- 
val < j < a; all the minus signs associated with the 
sublattice parity flips are removed from the spin correla- 
tions. 

Although the 'string' operator (— 1)^j ' is non- 
local it can be evaluated in a straightforward manner 
using the techniques introduced by Parola and SorcUa 
[Q. It is easily shown that in the large U limit, 

(M^(0) (-1)^^=1 ^""^ M'{x) ) 

x+l 

= Y.PsF{3){-^yOHe^sX3-l) (3) 

where Oueis. is the spin correlator of the Heisenberg 
chain, while P|^(j) = {n{0)n{x) 5{Y,'l^oni - j))sF is 
the probability of finding j spinless fermions in the in- 
terval [Q,x\. This factor causes the additional decay of 
the spin correlations due to the charge fluctuations in the 
standard spin correlator. However, it is easily shown that 
it is precisely compensated in Otop by the factor (— !)■' 
coming from the string operator and we find the result, 
asymptotically exact for large distances x. 
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Which is identical to the result for the pure Heisen- 
berg spin chain iln^^^(a;) after rescaling the amplitude 

of staggered spin M ^ and the measure of length 
X — > x/p, where p is the average charge density. Compar- 
ing Eq.(^ with the well-known 0| asymptotic behavior 
of the local staggered spin correlations for the Hubbard 
model, Oj^j{x) = (|Af^(0)M^(x)|) - """^j^'+T/^"^ we find 
that the former decays more slowly. The charge-stiffness 
Kp is associated with the decay of the charge correla- 
tions, {n{Q)n{x)) ^ cos([4fcF — ■n]x) / x'^'^ " , and Kp > 
for all p 7^ 1. Hence, the charge fiuctuations modify the 
spin correlations because the charge is attached to the 
sublattice parity flips, resulting in a simple multiplica- 
tive factor l/x^p . 

Is the above an accident of the strongly coupled case? 
Having identifled the operator Eq. (^ which in combi- 
nation with the local spin-spin correlator measures di- 
rectly the presence or absence of the topological order 
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it becomes possible to study it in any one dimensional 
system. For finite C/'s it becomes very hard to calculate 
Eq. from the Bethe-Ansatz solution and therefore 
we compute it numerically for the Hubbard model using 
the density matrix renormalization group. We employ 
an algorithm recently developed by one of us ||], mak- 
ing explicit use of the full non-abelian SO (A) symmetry 
of the Hubbard model. In this formulation, the gener- 
ators of the global symmetry group are the spin S and 
the pseudospin /; the latter is a generalization from the 
charge U{1) to an SU{2) symmetry such that the particle 
number at site i is given by = 1 -I- 21 f. 

In principle, it is straightforward to calculate the ex- 
pectation value of string operators such as Eq. (||) using 
the DMRG method. However, the calculation of criti- 
cal correlators with DMRG is not unproblematic even in 
1-l-lD, because eventually the truncation errors inherent 
in the form of the wavefunction will cause an exponential 
decay j^], and the requirement to use open boundaries 
causes a loss of translational invariance even at large dis- 
tances. We used a relatively large system size (1000 sites) 
to reduce the effect of the boundary conditions, and used 
a basis size (700 50(4) states) large enough to achieve 
a truncation length of the order of 200-300 lattice con- 
stants. In this way meaningful results can be obtained 
from a simple curve fit to obtain the desired exponents. 

In Fig. 2 we show our results for the exponents of both 
the standard two point spin correlator {M{0)M{x)), and 
our topological correlator, Eq. (|^), as a function of den- 
sity for various interaction strengths. The former illus- 
trates the quality of the calculation; the dominant corre- 
lation is at the wavevector 2kp — Tr, ~ cos{[2kp — tt]x) / x"^ 
and it is seen from Fig. 2 that the exponent rj depends 
strongly on the parameters. In fact, the exponent be- 
haves exactly according to the expectations: rj = K^+Kp 
where K„ = 1 and Kp is consistent with the values pre- 
viously obtained by Schulz [0 for the same values of 
U and p. Considering now the topological correlator, 
the dominant component lives at the wavevector q = 0; 
we find no other characteristic momenta in the Fourier 
transform of Eq. (^ which is already reminiscent of the 
staggered spin correlator of a Heisenberg chain. This is 
further amplified by our finding that the exponent rjtop 
does not depend on the microscopic parameters at all. In 
fact, for all parameters the exponent of the topological 
correlator rjtop = 1 = K^. Hence, regardless the values of 
U and p the long distance behavior of the topological cor- 
relator is indistinguishable from the spin-spin correlator 
of a Heisenberg chain, demonstrating that the Hubbard 
model indeed carries the sublattice parity topological or- 
der fully for all values U > 0. 

Let us now explain the reasons why bosonization fails 
in principle to describe the topological order. Tomonaga- 
Luttinger bosons cannot carry topological order for the 
elementary reason that these parametrize a relativistic 
theory. At the center of our considerations for the lat- 



tice model is the fact that the topological structure is 
Z2 valued (the sublattice parity) and this in turns rests 
on the fact that on a lattice the number operator Ui is 
integer valued: a site is either empty, singly or doubly 
occupied corresponding with charge quantum numbers 
0,1, and 2; by exponentiation one takes out the evenness 
or unevenness by the factors erf = (— 1)"' = ±1. In a 
relativistic theory one has to keep the fermion creation- 
and annihilation operator apart, — > ip'' {r)ip{r + e) and 
the resulting Schwinger terms are the building blocks for 
bosonization. Thereby, the charge density becomes the 
gradient of a continuous function n{x) — {1/Tr)dx- Ac- 
cordingly, (—1)"' ~ (—1)^=0 is a complex scalar and not 
an Ising valued operator. This reflects the fact it is im- 
possible to localize an electron in space in a relativistic 
theory, while the ability to ascribe a position to the elec- 
tron is at the heart of the squeezing procedure on the 
lattice. 
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FIG. 2. Exponents of the spin-spin correlation, rj, (top) 
and the topological correlation, rjtop, (bottom) as a function 
of electron filling, for various values of interaction U. 

In fact, this problem is to a certain extent repairable. 
The operator entering the string correlator is not it- 
self, but instead Hjcr'^ = exp(i7r dx n{x)), taking the 
continuum limit. The integral smooths out the discrete- 
ness, and an integral over the continuous bosonization 
fields can be manipulated to mimic the product of Ising 
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valued operators. This is a subtle UV regularization is- 
sue, which we will discuss elsewhere in detail. To give the 
reader some feeling, let us consider the expectation value 

of the string correlator itself, D{x) = ((— 1)^^=°"^) — 
{cos ttY, J n{j) ) -> {cos [tt dy n{y)]) . Accord- 
ing to bosonization, the total charge density is 
"(^) = \/f 1^0') + OcDwij) + O^cDwij) where 
OcDwix) — exp{~2ikFx) exp{i\/2Tr(pc) cos{\/2nips)- It 
follows that D{x) ^ {exp {^i\/2n [ipc{x) — fdO)])) — 
cos{2kpx)/x^p . On the other hand, because on the lat- 
tice (—1)"' = (— we might as well calculate D{x) — 

(cos 27r ^J^Q S'^(j) ) and using that in bosonization 

25^(2^) - + OsDwij) + oLvf(j) it follows 

that D{x) = For U ^ 0, ^ Kp and the 

two expressions are the same, but away from this point 
^ Kp and bosonizations runs into a paradox: de- 
pending the way one calculates D one obtains different, 
and mutually exclusive answers. In fact, the DMRG re- 
sult shows that D{x) = A/x^" + B cos{2kFx)/x^p {A, B 
are non- universal amplitudes), and apparently the two 
ways of calculating D{x) in bosonization recover part of 
the correct answer. A similar problem occurs when one 
uses bosonization to calculate the expectation value of 
nj(— l)^"' which is on the lattice just the identity op- 
erator; according to bosonization this decays like ^ittj! 
As we will discuss elsewhere ||ll|] , these problems can be 
traced back to spurious contributions associated with the 
'spreading' of the electron number density. 

In summary, we have discovered a correlation function 
which makes it possible to measure directly the presence 
of a hidden order underlying spin-charge separation in 
1-l-lD electron systems. This hidden order can be viewed 
as either a form of topological order, or as a geometri- 
cal structure where the spin system lives in a squeezed 
space different from the full chain. The suspicion was 
widespread that the Ogata-Shiba construction was spe- 
cial to the U oo limit and our main result is the numer- 
ical demonstration that squeezed space is generic in the 
scaling limit. We perceive this as a deep insight, powerful 
enough to hint at the existence of hitherto unidentified 
states of 1-l-lD matter. An exciting possibility is closely 
related to recent ideas regarding a potential connection 
between Ising gauge theory and the destruction of stripe 
order in 2+1 dimensions 12|. As we discussed, (—1)"* 
is Ising valued and the string operator can be written 
accordingly as ^icrf- this is just the Wilson line of Z2 
gauge theory Hence, if local Ising symmetry would 
be is dictating, {SUcr^S) ~ exp(a;/^g)(5'5') because the 
l.h.s. is gauge invariant while the r.h.s. is not; is the 
length scale where the gauge invariance emerges. In the 
Luttinger liquid the two are related by an algebraic factor 
1 /x^p but this is due to the algebraic order in the charge 
system. Therefore, under the conditions that (a) the spin 



system stays separated from the charge, (b) the charge 
quanta are bound to sublattice parity flips, and (c) the 
charge correlations become short range, one will obtain 
the required relations between the correlation functions 
signaling the Z2 gauge invariance. Condition (c) is gen- 
erally satisfied when the superconducting phase is forced 
to condense by applying an external Josephson field. The 
attractive U Hubbard model fails in this regard because 
its ground state can be viewed as a continuation of the 
local pair limit where the spin system is destroyed. How- 
ever, it appears |jl^ that the conditions might well be 
satisfied in t — J ladders under the influence of an exter- 
nal Josephson field, and these have in turn much more 
similarity with the situation in 2-l-lD. 

Acknowledgments. We acknowledge stimulating dis- 
cussions with S.A. Kivelson, S. Sachdev, F. Wilczek, T. 
Becker, E. Fradkin, N. Nagaosa and T.K. Ng. Financial 
support was provided by the Foundation of Fundamental 
Research on Matter (FOM), which is sponsored by the 
Netherlands Organization for the Advancement of Pure 
Research (NWO). Numerical calculations were performed 
at the National Facility of the Australian Partnership 
for Advanced Computation, through the Australian Na- 
tional University Supercomputer Time Allocation Com- 
mittee. 



[1] 
[2] 



[3] 
[4] 

[5] 
[6] 



[7] 
[8] 

[10: 
[11 

[12 



[13 
[14 



J. Voit, Rep. Prog. Phys. 58, 977 (1995) 
J. Zaanen O.Y. Osman, H.V. Kruis, Z. Nussinov and J. 
Tworzvdlo , Phil. Mag. B 81, 1 985 (2002); Z. Nussinov and 
J. Zaanen, ;ond-mat/0209437 ; J. Zaanen and Z. Nussinov, 



cond-mat/0209441 



E. H. Lieb and F.Y. Wu, Phys. Rev. Lett. 20, 1445 (1968). 

F. Woynarovich, J. Phys. C 15, 85 (1982); ibid. 15, 96 
(1982). 

M. Ogata and H. Shiba, Phys. Rev. B 41, 2326 (1990). 
This correlator was introduced first in the context of the 
dynamical stripes: J. Zaanen and W. van Saarloos, Physica 
C 282, 178 (1997). 

A. Parola and S. Sorella, Phys. Rev. Lett. 64, 1831 (1990). 
LP. McCulIoch and M. Gulacsi, Europhys. Lett. 57, 852 
(2002). 

S. Ostlund and S. Rommer, Phys. Rev. Lett. 75, 3537 
(1995). 

H. J. Schulz, Phys. Rev. Lett. 64, 2831 (1990). 

H.V. Kruis, LP. McCuIloch, Z. Nussinov and J. Zaanen, in 

preparation. 

Y. Zhang, E. Dernier and S. Sachdev, Phys. R ev. B 
66, 094501 ( 2002); S. Sachdev and T. Morinari, Font" 



mat/0207167 



J.B. Kogut, Rev. Mod. Phys. 51, 659 (1979). 

This appears to be the situation in doped t-J ladders, eg. 

S. R. White and D. J. Scalapino, Phys. Rev. Lett. 81, 3227 

(1998). 



4 



